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Abstract 


A  study  and  a  general  model  for  the  roll/pitch  frequency  "lock-in" 
phenomenon  are  made  for  antitank  long  kinetic  energy  projectiles 
and  missiles  in  general,  to  quantify  the  necessary  projectile 
damage/defect  required  to  cause  the  lock-in.  (Lock-in  is  a 
"persistent"  resonance  for  the  flight  duration  of  the  vehicle  [i.e., 
projectile].)  Fin  damage  and  body  mass  offset  are  modeled  as  the 
forcing  forces  and  moments  (asymmetries)  causing  the  lock-in. 
Idealized  fin  damage  is  modeled.  The  corresponding  pitch,  sideslip, 
and  roll  equations  are  numerically  integrated.  The  resulting  yaw  is 
largest  when  the  projectile  is  spinning  near  the  pitching  frequency, 
and  the  yaw  diminishes  when  the  spin  rate  (attributable  to  fin 
damage)  is  far  from  it.  The  nonlinear  induced  roll  moment  is 
essential  in  coupling  the  pitching  and  yawing  motions  with  the 
rolling  motion.  A  specific  and  complete  projectile  case  study  is 
presented,  in  which  quantified  fin  damage  and  rod  mass  offset  are 
found  to  cause  such  lock-in.  Computed  flight  motion  history  (pitch, 
side  slip,  and  roll),  and  yaw  amplification  factor  provide  insight  into 
this  lock-in  behavior.  The  model  is  general  andean  be  applied  to  any 
finned  projectile  or  missile. 
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FIN  DAMAGE  AND  ROD  ECCENTRICITY  FOR  SPIN/PITCH  LOCK-IN  FOR 
ANTIARMOR  KINETIC  ENERGY  PROJECTILES 


1.  INTRODUCTION 

Long  antitank  kinetic  energy  (KE)  projectiles  and  missiles  are  known  to 
occasionally  exhibit  abnormal  flight  behavior  characterized  by  (a)  spinning  at  a  different 
rate  than  their  intended  design  values,  and  (b)  having  larger  than  expected  yawing  motion 
angles.  In  fewer  cases,  this  behavior  can  be  catastrophic  to  the  motion  and  mission  of  the 
vehicle,  causing  an  ineffective  armor  penetration  for  the  KE  projectile  case,  or 
uncontrollable  flight  for  the  missile  case.  Note  that  the  projectile  flight  duration  is  only  1  to 
3  seconds. 

When  these  two  events  occur,  it  is  most  likely  that  a  roll/pitch  “lock-in”  has  occurred 
where  the  “vehicle”  (i.e.,  projectile/missile)  inappropriately  spins  because  of  damaged  fins  or 
other  asymmetry,  close  to  the  pitching  motion  frequency  of  the  vehicle.  This  fin  damage  is 
usually  referred  to  in  the  literature  as  “asymmetries,”  which  are  suffered  by  the  vehicle. 
Lock-in  is  a  “persistent”  resonance  for  the  flight  duration  of  the  vehicle.  Figure  1  shows  a 
fin  impact  imprint  indicating  damaged  parts  of  two  fins1,  bent  by  angles  exceeding  30°. 

Also,  a  KE  projectile  impact  imprint2  at  3  km  indicates  large  yaw  angle,  estimated  here  at 
about  8.8°. 

This  roll/pitch  lock-in  phenomenon  is  not  new  and  has  been  studied  by  many 
researchers  since  the  late  1940’s.  Nicolaides3  provided  an  early  analysis  while  Price4 
provided  a  useful  and  practical  model  for  the  mechanisms  causing  the  phenomenon. 

Clare6’6  provided  the  effect  of  nonlinearities  on  such  motion.  Murphy  included  the 
asymmetries  in  his  formulation1  and  later  showed8  that  lock-in  can  happen  at  a  reverse 
spin  rate.  More  recently,  Lin  et  al.9  studied  this  problem  for  a  re-entry  vehicle  (with  no 
fins)  and  considered  longitudinal  body  curvature  attributable  to  flexure  as  part  of  the  total 
angle  of  attack,  in  a  zeroth  order  model  for  flexure  effects.  For  more  prior  art,  the  reader  is 
referred  to  References  4,  5,  8,  and  9,  which  together,  list  about  186  references  covering  this 
lock-in  behavior  and  its  model  development  over  the  years. 

The  design  spin  for  a  KE  projectile  or  a  missile  has  to  avoid  being  close  to  the  natural 
pitching  motion  frequency  or  to  the  first  natural  lateral  bending  frequency  of  the  vehicle. 
For  the  KE  projectile  shown  in  Figure  2,  the  former  is  about  24  Hz  (decreases  from  25  to  23) 
while  the  latter  (for  the  steel  model  herein)  was  about  385  Hz.  The  design 
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SUBCALIBER  PROJECTILE 


DIMENSIONS  IN  CALIBER 


Figure  2.  Geometry  of  the  120-mm  M829  Subcaliber  KE  Projectile. 


spin  for  that  projectile  is  90  Hz,  as  depicted  in  Figure  3.  Certain  fin  damage  can  cause  the 
projectile  to  underspin  its  design  value  and  spin  near  its  pitching  frequency,  as  shown  in 
Figure  4.  These  particular  spin  values  near  this  pitching  motion  frequency  are  studied  in 
detail  later  in  this  report.  In  Reference  10,  the  spin  lock-in  with  the  upper  frequency  limit 
(the  first  bending  frequency)  attributable  to  structural  causes  was  studied  and  reported.  In 
the  present  lock-in  phenomenon,  the  mechanisms  for  causes  are  dynamic  motion  related  to 
the  pitching  and  side  slip  angles  coupled  to  a  lower  roll  rate  attributable  to  fin  damage. 

The  present  study  concerns  lock-in  at  this  lower  limit,  its  implication,  and  the  motion 
instability  amplification  it  induces.  Reference  11  provides  results  for  pitch  lock-in  with  rod 
flexure  but  gives  no  details  about  the  exact  formulation,  equations,  or  roll  information, 
including  any  induced  roll  moment. 

Usually,  “roll”  implies  a  low  spin  value.  Historically,  it  came  from  the  aeronautical 
reference  to  aircraft  roll  rate.  Spin,  however,  indicates  high  value  of  roll.  Commercial 
aircraft  may  roll  only  few  degrees  to  bank,  while  a  military  aircraft  may  do  a  few  complete 
rolls  for  maneuver.  Missiles  usually  have  low  roll  rate  (1  to  5  Hz),  while  finned  projectiles 
“spin”  at  50  to  150  Hz,  and  unfinned  projectiles  may  be  spun  to  more  than  300  Hz.  In  the 
present  work,  the  spin/pitch  lock-in  is  occasionally  called  roll/pitch  lock-in,  as  it  is 
frequently  called  in  the  literature. 

The  present  work  adapts  the  work  of  Price4  to  the  formulation  of  Reference  7. 
Numerical  solutions  will  be  obtained  for  the  three  coupled  roll/pitch/side  slip  equations.  A 
case  study  is  made  for  the  120-mm  M829  steel  version  of  the  projectile.  Flight  histories  and 
impact  imprints  are  presented  under  different  fin  damage  and  mass  offset  values. 

2.  ANALYSIS 

2.1  Definitions 

Resonance  between  the  roll  and  pitching  frequencies  is  defined  as  when  the 
former  is  affected  by  the  value  of  latter,  when  the  two  values  nearly  match  each  other.  Roll 
lock-in  is  when  the  former  frequency  attains  and  maintains  the  value  of  the  latter  for  a  long 
and  persistent  time  (1  to  3  seconds  for  KE  projectiles;  many  seconds  or  minutes  for  a 
missile).  Vehicles  can  roll  resonate  when  crossing  the  pitching  frequency  but  may  not  lock 
in  to  it.  Vehicles  that  roll  lock-in  must  resonate  with  the  pitch  frequency.  However, 
vehicles  may  temporarily  lock  in  but  then  escape  later,  in  what  might  be  called  “a  prolonged 
resonance.” 


4 


FREQUENCY  (Hz) 


120mm  M829  KE  PROJECTILE  (STEEL) 


400.0 


385 


300.0 


200.0 


LEGEND 

Pile  h  i  njg  .Motion  Fr  e  g  u  e  n  c  y . 

N a t u r g|  Ben d tng  Frequency _ 

ID  tended  Spin  De  sign  (Undomag  ed  F ins )_ 
Damaged  Fins,  36  Hz  S.S.  Spin 


100.0 


0.0 


Figure  3.  Design  Spin  (undamaged  fins)  Falling  Between  the  Two  Limiting  Frequencies. 
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Figure  4.  Spin  Cases  Studied  Near  the  Pitching  Frequency  (undamaged  fins). 


Yaw  angle  is  the  side  slip  angle,  J3.  The  yawing  motion  is  therefore  the  motion  in  the 
yaw  plane.  However,  it  is  common  in  ballistics  to  refer  to  the  total  angle  of  attack  as  the 
“yaw”  angle  and  the  total  motion  as  the  “yawing”  motion.  In  the  present  work,  yaw  refers 
to  the  correct  side  slip  definition. 

2.2  Induced  Roll  Moment 

The  induced  roll  moment,  C/j,  is  generated  on  a  finned  vehicle  because  of  the 
cyclic  fin  loading  at  different  angles  of  attack.  Its  magnitude  is  a  function  of  both  the  fin 
orientation  angle  during  roll  (with  respect  to  a  reference  plane)  and  with  the  total  angle  of 
attack.  It  is  nonlinear  (sinusoidal)  with  the  fin  orientation  angle  and  nonlinear  with  the 
total  angle  of  attack,  although  it  is  also  small  for  small  a.  More  importantly,  it  links  (i.e., 
couples)  the  yawing  motion  to  the  roll  motion.  Reference  4  only  cites  a  case  of  the  wind 
tunnel  experiment  of  the  Aerobee  150A  sounding  rocket,  but  the  actual  test  values  are 
provided  in  Reference  5  (which  provides  more  details  about  the  analytic  form  of  that 
coefficient).  C/j  is  defined  as 

C/j  =  induced  roll  moment/ (qArefd  sin(ny ))  (1) 

Clare5  provides  a  useful  algebraic  fit  based  on  the  experimental  testing  of  the  Basic 
Finner.  He  applied  the  same  expression  to  the  Apache  missile,  and  Price4  applied  it  to  the 
Aerobee  rocket. 

CH  =  (-0.87  +17.9  a*)at  (2) 

The  same  pattern  of  variation  with  at  was  also  observed  for  the  Basic  Finner  (code  name) 
projectile  and  is  given  in  Reference  12.  It  is  interesting  to  notice  that  both  cases  exhibit  a 
negative  value  until  at  =  12.6  and  then  become  rapidly  positive.  This  measured  behavior 
was  not  observed  or  used  in  Reference  8  which  uses  a  simplistic  C/j  algebraic  sixth  order 
variation  with  at  for  a  six-finned  projectile  configuration.  Clare  also  used  the  C/j  values  of 
the  Basic  Finner  for  his  Apache  missile  configuration,  which  is  generally  similar  in  shape  to 
both  the  Basic  Finner  and  the  M829  projectile  used  here.  In  the  present  work,  the 
expression  of  Equation  (1)  for  the  Basic  Finner  will  be  used  for  the  M829  projectile  since  no 
data  are  available  for  the  latter. 
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2.3  Mechanisms  for  Lock-in 


Reference  4  lists  three  main  missile  asymmetry  categories  that  can  cause  lock-in. 
First  is  configuration  asymmetries  such  as  fin  cant  manufacturing  tolerance,  fin  damage 
during  launch  or  in-flight,  or  a  bent  nose,  each  of  which  results  in  a  static  trim  angle. 
Rocket  thrust  misalignment  with  the  body  axis  is  another  example.  Second  is  the  mass 
center  location  deviation  from  the  body  axis,  which  provides  coupling  of  the  roll  with  the 
aerodynamic  normal  force  of  the  vehicle.  Third  is  the  aerodynamic- geometric  interactions 
such  as  the  fin  rotation-angle  of  attack  interaction  causing  the  induced  roll  moment 
described  earlier.  Elastic  bodies  that  can  bend  or  deform,  thus  causing  a  change  in  the 
aerodynamic  forces  from  their  intended  design  values,  are  another  example. 


2.4  Formulation 

2.4.1.  Analysis  of  Price  (Rolling  Coordinates  ) 

Price4  adapted  the  yaw  and  pitch  equations  developed  earlier  by  Nielson  for  a 
constant  spin  case,  to  accept  varying  spin  rates.  The  equations  of  Nielson  (cited  in 
Reference  4)  for  pitch  and  yaw  angles  were  written  in  body-fixed  coordinates,  i.e.,  rolling 
with  the  body  itself.  Price  then  added  the  roll  equation  and  included  the  induced  roll 
moment.  He  also  added  the  terms  for  the  effects  of  mass  offset,  moment  asymmetry,  and 
rocket  thrust  misalignment  to  the  three  equations.  The  latter  effect  will  not  be  considered 
in  the  present  work.  The  general  nomenclature  and  coordinates  are  shown  in  Figure  5. 

The  roll  motion  equation  with  the  added  mass  offset,  A  eg,  and  the  induced  roll  moment  was 
written  by  Price  as 

I X  =  q.Arefd.(Cjs.8  +  C\p  -  Cn.sin(n ~(j  )  -  .Qva-  a trim-sinQ -T ))  (3) 

V  d 


subject  to  the  initial  condition 


in  which 
( f=  d(f)  /  dt  =  p 

_l\  bp/cop  i 

Y  =  tan  - y - — — T  I 

_l.-p2  (l.  -lx/ fy)/ co 2p  J 


at  t  =  0, 


(j)  =  X  -  y 

b  —  (cjArej-  /  m  Vcop  ^CVf/  ji.  —  /Y  /  I f  J—  md  ~ C Mqf y 


(4) 
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Figure  5.  Nomenclature,  Coordinates,  and  Angles  (after  Price4). 


Also,  -  y/  is  the  phase  angle,  and  (f)  is  the  fin  orientation  angle.  The  sin( )  in  Equation  (3) 
represents  the  variation  with  the  roll  angle  for  each  of  Cu  and  the  mass  offset  roll  moment. 
Note  that  the  assymetry  pitching  moment  does  not  show  up  explicitly  in  the  roll  Equation 
(3),  except  through  the  Cu  term  which  reflects  the  values  of  a  and  |3,  which  in  turn  are 
affected  by  Cmo- 

2.4.2.  Non-Rolling  Coordinates 

The  analysis  of  Reference  7  is  almost  standard  in  ballistics  formulation.  It  uses  a 
ballistic  coordinate  system  (i.e.,  axes  are  attached  to  the  center  of  gravity  [CG]  of  the  body 
and  move  with  it  but  do  not  roll  with  the  body).  Reference  7  has  the  standard  form  which 
was  expanded  here  in  the  a,  j3,  and  time  variables,  modified  for  the  added  asymmetry 
effects,  and  re-written  in  Price’s  form: 


P 

a 


+  Aj  -  A2  a  -  Bj  J3  +  B2  a 
+  A]  a  +  A2  '(5  -  B\  a  -  B2  P 


=  Ci 

(5) 

=  c2 

(6) 

in  which 

A]  =  (qAref/ mV)(CLa  -  md2[CMq  +  Cy/d  ]/Iy) 

A2  =  -P(IX/Iy) 

B  l  =  (qAref  d/Iy)CMa 

B2  =  (pqA ref/m  V) ( C't a  VA’  +  md2/Iy  CMpa) 

C]  =  (qAref  d/I})(C Mo  cos  [A  +  tp]  -  C a  cos[r  +  </>]) 

C2  =  (qAref  d/Iy)(C Mo  sin  [A  +<j>]  -Ca^j  sin  [T  +(/)]) 

CMpa  =  Magnus  moment  slope  coefficient,  Magnus  moment/h/4 refd[pd/V] a) 

(CMq  +  Cm  a  )  —  pitch  moment  damping  coefficient,  pitch  damping  moment  l(qAref 
d[q  d/VJ) 

q  =  local  pitching  rate 

m  =  mass  of  projectile. 

The  K  angle  is  such  that  Cmoz  ~  Cmocos0-)  and  C,y/or  =  C,y/0.s'hi(A) .  The  pitch  angle,  a,  is 
positive  when  the  vehicle  nose  is  up,  while  Cm  (hence,  CmP  is  negative  when  it  turns  the 
nose  tip  down.  Without  C;  and  C2  in  Equations  (5)  and  (6),  both  a  and  p  equations  will 
have  cop  as  their  natural  frequency. 
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The  roll  equation  [Equation  (3)]  is  modified  for  the  non- spinning  axes  and  re-written 
after  Price’s  model.  The  modification  is  to  reflect  the  tracking  of  the  roll  angle  in  non¬ 
rolling  coordinates.  The  angles  of  rotation  are  tracked  from  a  reference  plane  taken  to  be 
the  vertical. 

Ix  t  =  q-Aref.d.(Cfo.b  +  Cjp  -  Cn.sin(ruj) )  -  .Ch/U.  c it.sin[<f>  +  T -(f)  ]) 

V  d 

(7) 

in  which 

!>  =  (</>+  <t>0  -  <t>  ) 

(f)  —  roll  angle,  measured  from  a  fixed  reference  plane 

(f>0  =  initial  roll  angle,  measured  from  same  reference  plane 

A 

(f)  =  tan^(j3/  a),  pitching  motion  plane,  measured  from  same  reference  plane 


in  which  (f>  is  the  the  angle  between  the  pitch  plane  and  the  fin  orientation  angle.  The  yaw 
amplification  factor  for  the  trim  angle  of  attack,  based  on  Reference  7  for  symmetric 
missiles  with  small  yaw  angles  and  linear  aerodynamics,  is  usually  written  as13 

—0.5/?  Arej-d  Cfy[a 


AMP  =  ■ 


a, 


a 


trim  at  zero  spin 


({p  -  ¥f  }{p  -¥s  }-^fAs^  +  (As  {/?  -  (ff  }+  Af  {/, 


in  which 


p 

=  p/V 

=  plx(l.  +  a)/(2Iy) 

*s 

=  plx(l.  -  a)/ (21  y) 

a 

=  Jl.-1./Sg 

Sg 

t2— 2 

-  2i,p 

/  (nIypCMad3  ) 

Af 

II 

[-CjVa(l.-l./CT)+0.5*2l 

^1. + 1.  /  crJpMq  +  kj  CMpa  /  crj 

As 

*  * 
II 

[-CKo(l.+  l./cr)+a5^l 

(l . — 1 . / (7^jZMq  —  k/  C'\fpa  /  crj 

hi 

=  md2/Ix 

h2 

=  md2/Iy 
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This  amplification  factor  represents  the  ratio  of  the  steady  state  yaw  angle  of  a  spinning 
projectile  to  the  corresponding  steady  state  yaw  angle  of  a  corresponding  non-spinning 
projectile  (i.e.,  projectile  in  combined  pitching  and  yawing  motions  only).  In  the  present 
study,  however,  the  time-dependent  growth  history  in  the  yaw  angle  (in  reference  to  the 
initial  total  yaw  angle  at  launch)  is  introduced  in  Section  5.  This  new  amplification  factor 
time  history  is  also  given  in  Section  5. 

2.4.3  Fin  Damage 

In  ballistics  practice,  damage  to  fins  may  be  observed  on  x-ray  or  video  films  taken 
for  projectiles  exiting  the  muzzle.  Some  fin  damage  is  only  realized  when  yaw  tests  (using 
yaw  cards)  are  performed.  However,  fin  damage  cannot  be  predicted  a  priori.  Most  fin 
damage  occurs  in  the  gun  tube  and  before  the  projectile  exits  the  muzzle  (some  damage 
because  of  aeroheating  may  occur  during  flight),  as  was  also  recently  ascertained  in 
Reference  14.  Seldom  are  the  fins  recovered  after  each  firing  to  examine  the  damage.  In 
many  cases,  fin  damage  could  only  be  assumed,  based  on  untypical  flight  performance  such 
as  large  yaw  angles  at  large  ranges  (1  to  3  km  for  projectile  applications)  or  a  screen  impact 
imprint  indicating  rod  bending  or  even  missing  or  bent  fins.  Unfortunately,  for  those 
observed  flight  abnormalities,  no  fin  damage  was  expected  nor  were  the  fins  recovered. 
Therefore,  theoretical  studies  and  models,  like  the  present  work,  are  used  to  explain  the 
relation  between  the  abnormal  flight  behavior  (e.g.,  roll/pitch  lock-in  or  roll/flexure  lock-in) 
and  possible  fin  damage.  A  question  is  often  posed  about  how  much  and  what  type  of  fin 
damage  it  takes  to  cause  such  behavior,  even  if  the  damage  itself  cannot  be  prevented.  The 
present  work  sheds  some  light  in  that  respect.  Theoreticians8  assign  numerical  values  to 
“asymmetries,”  but  they  never  relate  these  to  real-life  fin  damage  and  examine  if  those 
values  are  realistic. 

In  the  present  study,  fin  damage  is  idealized  into  three  types  that  are  most  often 
observed  from  yaw  card  tests.  Corresponding  damaged  fin  areas  of  5%  and  10%  for  a  single 
fin  panel  bent  at  different  damage  angles  are  studied  as  a  practical  example.  The  fin 
damage  cannot  be  just  arbitrarily  assumed  in  size  and  direction  but  must  be  made 
consistent  to  explain  the  loss  in  roll  rate  from  its  design  value  as  well.  Therefore,  only 
damage  in  a  certain  direction  relative  to  the  fin  cant  angle  and  spin  direction  can  be 
applied.  Although  other  fin  angle  damage  directions  obviously  happen,  they  will  not  result 
in  lower  roll  and,  thus,  in  this  spin/pitch  lock-in.  In  fact,  Reference  10  purposefully 
considered  the  reverse  case  only  (i.e.,  only  damage  that  will  overspin  the  projectile)  to  study 
the  spin/flexure  lock-in  behavior.  Figure  6  shows  the  required  fin  damage  direction  to 
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Figure  6.  Leading  and  Trailing  Edge  Idealized  Fin  Damage. 


o-> 


cause  a  change  in  fin  normal  force  and  roll  moment  for  an  equivalent  fin  panel  canted  at 
8 eg.  The  projectile  roll  is  positive  when  clockwise  (looking  forward  from  the  tail),  and  8eg  is 

positive  when  it  is  up  in  the  direction  of  producing  positive  roll.  The  static  (non-rolling) 
asymmetries  A  Cjy  and  A  Cm  (Cmo)  can  be  positive  or  negative,  depending  on  the  damage 
location  and  direction.  Further  details  are  discussed  in  the  case  study  section. 

2.4.4  Mass  Offset 

The  projectile  radial  mass  asymmetry  is  represented  by  the  mass  offset  radial 
distance,  A  eg.  This  offset,  which  always  exists  in  real  life  with  different  degrees  of  severity, 
creates  a  rolling  moment  because  of  the  normal  force  that  acts  on  the  aerodynamic,  not  the 
mass-center  point  of  the  rod  cross  section.  This  offset  also  gives  rise  to  pitching  and  yawing 
moments  attributable  to  the  aerodynamic  axial  force  (C^),  as  can  be  seen  in  Figure  5  and 

reflected  in  Equations  (5)  and  (6).  An  axial  mass  shift  (i.e.,  shift  along  the  body  axis)  was 
not  considered  explicitly  as  an  asymmetry  moment  in  this  work.  However,  it  can  be 
considered  as  part  of  the  off-design  asymmetry  moment,  M0,  used  here  and  which  has 

primarily  focused  on  moments  resulting  from  fin  damage. 

A  value  of  A  eg  =  0.20d  is  extremely  large  and  has  been  used  only  as  an  application 
example.  Actual  practical  values  are  of  the  order  of  0.003d  to  0.009d  for  projectile 
applications.  For  missile  applications,  this  offset  is  usually  larger  than  those  given  for 
projectiles. 

3.  CASE  STUDY 

Calculations  were  performed  to  simulate  the  case  of  a  subcaliber  projectile  that  has  a 
length/ diameter  (L/d)  of  22.72,  as  given  in  Figure  2.  The  aerodynamic  coefficients  were  obtained 
from  the  U.S.  Army  Research  Laboratory  (ARL)  enclosed  transonic  range  test  database.  In  the 
enclosed  range,  the  tests  are  made  on  steel  models  rather  than  on  the  actual  design  rounds  for 
safety  considerations.  Therefore,  simulations  were  also  made  for  that  particular  case.  The 
relevant  information  regarding  this  configuration  is  given  in  Table  1. 


10 


Table  1.  Physical  Properties  of  the  Projectile 


Ref.  diameter 

— 

1.06  in. 

(25.05  mm) 

L/d 

= 

22.72 

-- 

Material 

= 

steel 

-- 

Mass 

= 

4.281  lb 

(1.946  kg) 

lx 

= 

0.6956  lb-in2 

(2.04xl0"4  kg-m2) 

h 

= 

143.56  lb-in2 

(0.42xl0_1  kg-m2) 

Launch  speed 

= 

5905.8  ft/s 

(1806  m/s) 

Launch  spin 

— 

0  Hz 

"" 

Twenty-three  cases  with  four  different  steady  state  spin  values  were  computed,  covering 
the  spectra  of  above/at/below  the  pitching  frequency  case.  The  pitching  frequency  was  25  Hz 
at  launch  (M  =  5)  and  dropped  to  23  Hz  at  Mach  =  4.0  (or  about  t  =  0.9  s).  The  different 
steady  state  spin  cases  were  considered  to  reflect  the  effect  of  damaged  fins.  The  actual 
equivalent  fin  cant  angle  for  the  projectile  is  estimated15  to  be  0.55°,  resulting  in  a  steady 
state  spin  of  90  Hz.  Table  2  provides  the  fin  cant  angle  values  applied  to  simulate  damaged 
fins  and  the  corresponding  steady  state  spin  values.  These  chosen  spin  values  were  selected 
so  as  to  provide  insight  about  resonance  behavior  near  and  through  the  pitching  frequency, 

(Op. 


Table  2.  Equivalent  Fin  Cant  and  Steady  State  Spin 


Case 

No. 

Fin 

Status 

Equiv.  Fin 
Cant  Angle, 
(deg) 

Cant  Angle 
as  Fraction 
of  5d 

Steady  State 
Spin  (Hz) 

I 

undamaged 

0.55 

(1.0  Sd) 

90 

Ha 

damaged 

0.22 

(0.4  Sd) 

36 

lib 

damaged 
(spin  above  cop) 

0.165 

(0.3  Sd) 

27 

HH 

1 — 1 

1 — 1 

damaged 
(spin  at  cop) 

0.137 

(0.25  Sd) 

22.5 

IV 

damaged 
(spin  below  (Op) 

0.11 

(0.2  Sd) 

18 
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For  quantified  fin  damage,  three  idealized  types  were  modeled,  reflecting  leading, 
trailing,  and  tip  edge  damage,  as  shown  in  Figure  7  for  the  six-fin  configuration.  The  actual 
fin  shown  in  Figure  8  is  modeled  by  an  equivalent  fin  panel  wholly  canted  at  8eq  =  0.55°. 

Representative  fin  damage  surface  areas  for  the  three  types  are  also  shown  in  the  same 
figure.  The  corresponding  change  in  fin  normal  force,  A  Cjy,  and  A  Cm  were  computed  and  are 

presented  in  Table  3.  An  idealized  equivalent  uniform  fin  cant  angle  damage  for  the  six  fins  is 
listed.  A  single  panel  with  a  damaged  area  equivalent  bend  angle  is  modeled  to  provide  the 
same  lower  spin  rate  that  the  six  damaged  fins  would  have.  The  corresponding  change  in  the 
fin  normal  force,  A Cjy,  and  the  corresponding  change  in  Cm  (be.,  Cmo )  were  then  computed. 
The  corresponding  values  for  damage  Types  I,  II,  and  III  are  all  listed  in  Table  3.  The  C]ya  of 
the  fin  set,  including  its  body  interference  effect,  was  computed  by  using  the  Naval  Surface 
Warfare  Center  aerodynamic  prediction  (NSWC-AP95)  fast  design  code16-  Linear 
aerodynamics  were  used  to  provide  engineering  estimates  for  the  fin  damage  angles  in 
association  with  the  damaged  fin  panel  areas.  For  example,  a  5%  area  fin  damage  bent  at  8° 
would  be  equivalent  to  a  10%  fin  area  damage  bent  at  40°.  This  is  based  on  the  relation  that  a 
fin  panel  lift  is  linearly  proportional  to  the  product  of  its  surface  area  and  the  surface  angle  of 
attack.  Although  this  linearity  may  not  be  valid  for  large  angles,  it  simplifies  the  procedure. 
More  sophisticated  nonlinear  models  may  obviously  be  devised. 
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Figure  7.  Three  Types  of  Idealized  Fin  Damage. 
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Figure  8.  Quantified  Fin  Area  Damage  for  Types  I,  II,  and  III. 


Table  3.  Damaged  Fin  Areas  and  Angles  With  Their 
Corresponding  A  Cat  and  A  Cm 


Resulting 

Reduced 

Roll  Rate 

(Hz) 

Uniform  Damage 

to  Six  Fins 

Equivalent  Damage 

to  a  Single  Fin 

Single  Fin  Damage  Parameters  and  Forces 

Damage 

Type  I 

Damage 

Type  n 

Damage 

Type  Ed 

Initial 

(deg) 

5d, 

(deg) 

Final 

6„ 

(deg) 

j  |  i  ! 

Damaged 

Area 

(Ad/Af) 

8* 

(deg) 

A1 

(deg) 

51 

(deg) 

A2 

(deg) 

52 

(deg) 

A3 

(in) 

acM  (^Mo) 

Damage  Type 

1 

D 

m 

36 

0.55 

-0.33 

0.22 

-1.43 

0.05 

-28.6 

16 

+28.6 

1.5 

-28.6 

28.6 

-0.094 

-0.97 

+0.78 

+0.93 

0.10 

-14.3 

30 

+14.3 

3.0 

14.3 

27 

0.55 

-0.385 

0.165 

-1.76 

0.05 

-35.2 

16 

+35.2 

1.5 

-35.2 

0.(J7h 

-35.2 

-0.116 

+  1.20 

tO.96 

+1.15 

0.10 

-17.6 

30 

+17.6 

3.0 

-17.6 

0.1 6h 

-17.6 

22.5 

0.55 

-0.412 

0.138 

-1.92 

0.05 

-38.4 

16 

+38.4 

1.5 

-38.4 

0.07h 

-38.4 

-0.126 

+  1.31 

+1.04 

+1.25 

0.10 

-19.2 

30 

+19.2 

3.0 

-19.2 

18 

0.55 

-0.44 

0.11 

-2.09 

0.05 

-41.8 

16 

+41.8 

1.5 

-41.8 

0.07h 

-41.8 

-0.137 

+1.42 

+1.14 

+1.36 

0.10 

-20.9 

3.0 

-20.9 

0.1 6h 

-20.9 

Different  independent  values  for  the  mass  offset  were  used,  ranging  between  O.Od  and 
0.20d.  C^a,  for  the  undamaged  projectile  was  -35  at  the  launch  speed.  Values  for  Cmo 
varying  between  0.0  and  +2.444  (corresponding  to  0.0°  and  -4.0°  static  trim  angle)  were  used. 
The  numerical  values  used  for  both  parameters  are  listed  separately  in  Table  4. 
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Table  4.  List  for  the  Independent  Values  Used  for  Mass 
Offset  and  Asymmetry  Moment 


Mass  Offset, 
Acg/d 

CMo 

The  Corresponding  Trim 
Angle  for  Cmo  (deg) 

0.0 

0.0 

0.0 

0.001 

+0.061 

-0.1 

0.004 

+0.122 

-0.2 

0.01 

+0.305 

-0.5 

0.04 

+0.428 

-0.7 

0.08 

+0.489 

-0.8 

0.10 

+0.916 

-1.5 

0.12 

+1.222 

-2.0 

0.20 

+2.444 

-4.0 

4.  COMPUTATIONS 

Computations  were  made  by  solving  the  three  equations  as  an  initial  value  problem. 

The  alpha  and  beta  equations  were  solved  simultaneously,  then  followed  by  the  roll  equation. 
Small  time  steps  were  used  (At  =  3.0xl0'6  s).  Three-point,  central  differences,  second  order 
accurate  (in  time)  finite  differencing  was  used,  with  recurrence  relations  relating  the  variable 
value  at  the  n+1  time  level  to  the  previous  n  and  n-1  time  levels.  Knowing  the  first  two 
values  at  time  zero  and  at  At  using  the  initial  conditions,  one  can  then  proceed  directly. 

The  roll  equation  was  solved  as  a  first  order  equation  in  ,  with  the  boundary 
condition  given  by  Equation  (4).  Equation  (2)  was  written  as 

+  0„+Ml»  £+cd„=0  (8) 

The  solution  was  then  computed  as 

¥,*l  =  &  -  Ar  [ss|  „  t,  +  cc|„  (9) 

The  roll  equation  is  coupled  to  the  pitch/yaw  equations  through  the  induced  roll 
moment  term  which  includes  at  and  C/p 

Calculations  were  usually  stopped  at  about  t  =  0.9  s,  which  corresponds,  for  the  case 
considered,  to  a  distance  of  about  5,578  ft  (1700  m)  down  from  the  gun  muzzle.  This 
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corresponds  to  about  300,000  time  steps.  A  single  case  takes  about  15  seconds  on  a  UNIX™ 
minicomputer  (workstation). 

5.  RESULTS 

Because  of  the  many  parameters  of  the  problem,  one  set  of  similar  conditions  was 
applied  so  that  the  spin/pitch  lock-in  phenomenon  can  be  observed  without  changing  the 
parameters.  Later  studies  may  be  made  with  these  parameters  as  variables  to  discuss  their 
individual  influences  in  detail. 

The  initial  conditions  were,  for  example,  held  the  same  for  all  cases.  The  values  a(t0) 

=  4°;  /? (t0)  =  3°;  <^t0 )  =  1.0  rad/s;  and  >^(t0)  =1.0  rad/s  were  used.17  The  roll  equation  was 
also  solved  under  zero  initial  values  for  both  (f)  and  .  Any  other  values  may  obviously  be 

prescribed  for  other  desired  conditions.  Relatively  large  values  for  the  initial  yaw  angles  (a 
=  4°  and  p  =  3°)  were  used  to  illustrate  the  point  that  these  angles  will  quickly  dampen  to 
small  values  close  to  zero  in  no-lock-in  cases.  This  usually  occurs  after  less  than  0.6  second, 
or  about  800  meters  from  the  the  muzzle,  for  the  present  configuration. 

Other  parameters  are  the  mass  offset  initial  location  angle,  T ,  and  the  Cmo  moment 
orientation  angle,  A  .  They  were  taken  to  be  30°  and  70°,  respectively.  Reference  4 
provides  a  discussion  about  A  angle  role  in  the  lock-in.  A  list  of  cases  computed  in  this 
study,  including  the  corresponding  values  of  the  parameters  used,  is  provided  in  the  form  of 
Table  5. 

First,  the  case  of  spinning  above  cop  with  steady  state  spin  of  36  Hz  is  presented.  This 
case  Number  1,  with  no  fin  damage  or  asymmetries  was  computed  and  the  typical  a  -  and  /?- 
angle  histories  are  given  in  Figure  9,  both  approaching  zero  values  at  about  t  =  0.6  s. 

Case  Numbers  2  through  10  were  computed  with  varying  A  eg  but  with  zero  Cmo ■  For 
the  higher  spin  rate  of  36  and  27  Hz,  no  lock-in  occurred  with  the  given  values  of  A  eg,  as 
shown  in  Figures  10  and  11,  respectively.  Although  the  value  0.20d  is  relatively  large  and 
is  not  expected  in  practice,  one  must  remember  that  no  Cmo  was  applied  here,  which  is  not 
the  case  in  real  life.  However,  as  the  steady  state  spin  caused  by  damage  gets  closer  to  the 
pitch  frequency,  such  as  for  the  cases  of  22.5  and  18  Hz,  lock-in  occurs  with  a  smaller  value 
of  A  eg.  Figure  12  indicates  the  lock-in  for  the  18-Hz  case  after  about  0.5  s  from  launch. 
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Figure  9.  Pitch  and  Side  Slip  Angles  for  No-Lock-in  Case 


FREQUENCY  (Hz) 


Figure  10.  Roll/pitch  Resonance  and  Then  Escaping  Pitch  Lock-in  ( 8eq  =  0.4  8^  ,  Cmo  =  0.0,  and  increasing  A  eg). 


FREQUENCY  (Hz) 


Figure  11.  Roll/pitch  Resonance  and  Then  Escaping  Pitch  Lock-in  ( 8eq  =  0.3  8^  ,  Cmo  =  0.0,  and  increasing  A  eg ). 


FREQUENCY  (Hz) 


Figure  12.  Roll/pitch  Resonance  and  Roll/pitch  Lock-in  (8e„  =  0.2  8^  ,  Cmo  =  0.0,  and  increasing  A  eg ). 


Table  5.  Cases  Computed  and  Parameter  Values 


Case 

No. 

Equiv.  Fin 
Cant  Angle 

Cmo 

Cmo  Correspond. 
Trim  Angle  (deg) 

A  cg/d 

Lock-in 

Status 

1 

0-4  §d 

0.0 

0.0 

0.04 

no  lock 

2 

0-4  8d 

0.0 

0.0 

0.08 

no  lock 

3 

0-4  8d 

0.0 

0.0 

0.20 

no  lock 

4 

0-3  8d 

0.0 

0.0 

0.0 

no  lock 

5 

0-3  Sd 

0.0 

0.0 

0.04 

no  lock 

6 

0-3  8d 

0.0 

0.0 

0.12 

no  lock 

7 

0.25  5d 

0.0 

0.0 

0.0 

no  lock 

8 

0.25  8d 

0.0 

0.0 

0.08 

lock-in 

9 

0-2  Sd 

0.0 

0.0 

0.0 

no  lock 

10 

0-2  Sd 

0.0 

0.0 

0.10 

lock-in 

11 

0-4  8d 

+0.305 

-0.5 

0.0 

lock-in 

12 

0-4  8d 

+0.916 

-1.5 

0.0 

lock-in 

13 

0-4  8 d 

+2.444 

-4.0 

0.0 

lock-in 

14 

0-3  8d 

0.0 

0.0 

0.0 

no  lock 

15 

0-3  Sd 

+0.428 

-0.7 

0.0 

lock-in 

16 

0-3  Sd 

+1.222 

-2.0 

0.0 

lock-in 

17 

0.25  5d 

0.0 

0.0 

0.0 

no  lock 

18 

0-25  Sd 

+0.061 

-0.2 

0.0 

lock-in 

19 

0-2  Sd 

0.0 

0.0 

0.0 

no  lock 

20 

0-2  8d 

+0.610 

-2.0 

0.0 

lock-in 

21 

0.4^ 

+0.061 

-0.10 

0.001 

no  lock 

22 

0-4  <5^ 

+0.122 

-0.20 

0.040 

no  lock 

23 

0-4  Sd 

+0.489 

-0.80 

0.100 

lock-in 

Case  Numbers  11  through  20  examine  the  effect  of  increasing  the  value  of  Cmo 
without  the  mass  offset  influence  being  present.  Figure  13  depicts  the  effect  of  increasing 
Cmo  for  the  fin  damage  spin  of  36  Hz.  It  is  interesting  to  note  that  for  smaller  fin  damage, 
it  takes  longer  to  lock  in,  while  for  larger  damage,  lock-in  occurs  faster  and  with  closer 
value  to  the  pitch  frequency.  The  same  behavior  is  depicted  in  Figure  14  for  the  fin  damage 
of  27-Hz  steady  state  spin.  Figure  15  depicts  the  yaw  angle  amplification  factor  for  this 
later  case.  The  amplification  factor  here  is  defined  as  the  computed  total  angle  of  attack, 
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Figure  13.  Roll/pitch  Lock-in  ( 8eq  =  0.48^,  increasing  Cmo  >  Acg  =  0.0). 


FREQUENCY  (Hz) 


Figure  14.  Roll/pitch  Lock-in  and  No  Lock-in  ( 8eq  =  0.4  8^  ,  increasing  Cmo  >  Acg  =  0.0). 


YAW  AMPLIFICATION  FACTOR 


120mm  M829  KE  PROJECTILE  (STEEL) 


Figure  15.  Yaw  Angle  Amplification  for  Roll/pitch  Lock-in  ( 8eq  =  0.3  8f/ ,  increasing  Cmo  >  Acg  =  0.0). 


divided  by  the  initial  total  angle  of  attack.  Amplification  factors  of  1.2  and  3.6  are  noticed 
for  cases  of  lock-in  with  fin  damage,  while  attenuation  to  a  value  of  zero  is  noticed,  as 
expected,  for  the  case  with  no  fin  damage  (no  lock-in). 

In  real  life,  usually  cases  with  both  asymmetries  and  mass  offset  occur.  Cases  21 
through  23  were  computed  showing  increased  values  for  both  Cmo  and  A  eg.  Pitch 
resonance  changes  to  pitch  lock-in  when  both  parameters  were  increased  as  reflected  by 
Case  Numbers  21  through  23  of  Table  5  and  shown  in  Figure  16.  Figure  16  also  shows  the 
change  from  short-lived  pitch  resonance  to  full  lock-in  with  cop  when  the  asymmetry 
damage  increased.  Figure  17  shows  the  a  -  and  (3 -angle  histories  for  Case  Number  21, 
which  reflects  very  small  asymmetries.  The  plot  indicates  almost  unaffected  motion. 

The  expected  impact  imprint  on  cardboard  yaw  cards  (or  wire  mesh  screens)  used  in 
measuring  the  “yaw”  motion  in  practical  ballistics  can  now  be  constructed  from  the  results 
computed.  Figure  18  presents  the  a  -  and  (3  -angle  histories  for  Case  Number  12  for  lock-in. 

Figure  19  depicts  the  expected  imprint  for  the  lock-in  case  of  Figure  18  at  t  =  0.633  s,  in 
which  |a|  and  |/?|  were  computed  to  be  3.54°  and  1.28°,  respectively,  with  a  total  yaw  angle 

of  3.76°.  A  harsher  case  of  lock-in  with  a  larger  Cmo  is  reflected  in  Figures  20  and  21.  The 
| O']  and  |/?|  computed  were  9.44°  and  5.60°,  respectively,  thus  yielding  a  total  yaw  angle  of 

11.97°  at  t  =  0.194  s. 

In  the  practice  of  ballistics,  the  a  -  f3  plot  is  usually  monitored  to  provide  insight  into 
the  flight.  The  three  cases  of  no  lock-in,  late  lock-in,  and  immediate  lock-in  are  presented 
in  Figures  22  through  24.  Figure  22  provides  the  typical  KE  flight  pattern  with  no 
asymmetry  and  therefore,  no  lock-in.  Figure  23  shows  Case  Number  11  which  locked  in 
later  in  its  flight,  as  was  also  indicated  by  Curve  Number  1  of  Figure  13.  Figure  24  depicts 
the  case  of  a  quick  roll/pitch  lock-in  of  Case  Number  23.  All  cases  started,  as  can  be  seen, 
with  the  same  initial  conditions  of  a  =  3°  and  /?  =  4°. 

Additional  cases  can  now  be  routinely  computed.  Reference  4  provides  a  simple 
“guideline”  for  the  expected  values  for  both  asymmetry  parameters  for  lock-in.  These 
values  may  be  used  as  first  estimate  for  the  value  of  the  parameter  under  which  the  vehicle 
may  lock  in.  The  exact  value  can  then  be  determined  by  computing  the  case  with 
perturbations  on  these  estimates  and  monitoring  the  resulting  spin  history,  and  the  a  —  [3 

plot. 
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Figure  16.  Roll  Resonance  and  Roll/pitch  Lock-in  ( 8eq  =  0.4  8^  and  increasing  Cmo  and  Acg ). 
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Figure  17.  Pitch  and  Side  Slip  Angles  for  a  Case  of  Very  Slight  Asymmetries  ( 8eq  =  0.4  8^,  Cmo  =  +0.061,  and  A  eg  =  0.001). 


PITCH  and  SIDE  SLIP  Angles  (deg) 


Figure  18.  Pitch  and  Side  Slip  Angles  for  a  Case  of  Roll/pitch  Lock-in  ( 8eq  =  0.4  8^  ,  Cmo  =  +2.444,  and  A  eg  =  0.0). 


Figure  19.  Predicted  Impact  Imprint  of  a  Case  of  Roll/pitch  Lock-in  (8eq  =  0.4  8^  ,  Cmo  =  +0.916,  and  A  eg  =  0.0). 


PITCH  and  SIDE  SLIP  Angles  (deg) 


Figure  20.  Pitch  and  Side  Slip  Angles  for  a  Case  of  Roll/pitch  Lock-in  ( 8eq  =  0.3  8^  ,  Cmo  =  +1-222,  and  A  eg  =  0.0). 


Figure  21.  Predicted  Impact  Imprint  of  a  Case  of  Roll/pitch  Lock-in  ( 8eq  =  0.3  8^  ,  Cmo  =  +1-222,  and  A  eg  =  0.0). 


Figure  22.  a  -  J3  Plot  for  a  No-Lock-in  Case  ( 5eq  =  0.4  8^  ,  Cmo  =  -0.0,  and  A  eg  =  0.0). 


ALPHA  (deg) 


Figure  23.  a  -  fl  Plot  for  a  Case  With  Roll  Resonance  and  Later  Roll/pitch  Lock-in  ( 8eq  =  0.4  8^  ,  Cmo  =  +0.305,  and  A  eg  =  0.0). 


Figure  24.  a  -  (3  Plot  for  a  Lock-in  Case  ( 5eq  =  0.4  8^  ,  Cmo  =  +0.489,  and  A  eg  =  0.10). 


6.  CONCLUSIONS  AND  SUMMARY 


The  following  conclusions  are  supported  by  the  results  obtained  through  the 
application  of  the  present  analysis  and  model: 

1.  Fin  area  damage,  in  certain  size  and  direction,  can  cause  the  lowering  of  the  design 
steady  state  spin  to  a  value  close  to  the  vehicle  pitching  motion  frequency,  thus  creating  the 
environment  for  possible  locking  in  to  that  frequency. 

2.  It  does  not  take  a  large  fin  damage  to  cause  the  low  spin  and  therefore  provides  the 
condition  for  possible  locking  in  at  the  pitching  motion  frequency.  For  the  case  studied  for  a 
six-finned  projectile,  a  single  fin  panel  idealized  area  damage  of  10%  of  the  fin  area,  bent  at 
19.2°,  was  enough  to  lower  the  spin  and  also  cause  the  lock-in. 

3.  The  higher  or  lower  the  expected  steady  state  spin  (attributable  to  damage)  differs 
from  the  pitching  motion  frequency,  the  larger  the  mass  offset  needed  to  cause  lock-in,  in 
absence  of  fin  asymmetry  damage. 

4.  If  the  fin  damage  is  less  than  that  required  to  cause  lock-in,  the  spinning  motion 
will  escape  from  being  locked  in  but  with  some  resonance  effect  on  the  spin  history. 

5.  In  case  of  mass  offset  alone,  lock-in  can  also  theoretically  happen  with  no  moment 
asymmetries  (Cmo  =  0.0). 

6.  If  a  combination  of  both  asymmetry  moment  and  mass  offset  is  causing  the  lock-in, 
their  individual  required  magnitudes  will  be  less  than  those  required  if  each  were  to  have 
been  applied  alone  to  cause  a  lock-in  (confirming  one’s  expectation). 

In  summary,  a  model  was  developed  to  study  and  simulate  the  roll  and  pitching 
motions  lock-in  observed  for  some  projectiles  in  flight.  Fin  damage  and  mass  offset  were 
quantified  and  their  effects  on  the  motion  of  the  vehicle  were  computed.  Flight  histories  of 
the  motion  including  the  spin  rate,  pitch  and  yaw  angles,  and  yaw  amplification  factor  are 
computed.  Both  roll  resonance  and  roll  lock-in  behavior  were  reflected  in  the  flight 
histories  computed.  Projectile  impact  imprint  on  yaw  cards  is  predicted  and  simulated  for 
any  distance  from  the  muzzle.  The  model  established  is  applicable  for  missiles  or  finned 
projectiles.  A  case  study  was  performed  in  detail  for  a  specific  projectile,  as  a  real-life 
applied  example. 
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LIST  OF  SYMBOLS 


Ad  =  damaged  area  of  a  fin  panel,  in.2 

Af  =  area  for  one  single  fin  panel,  in.2 

Aref  =  reference  area  (^"d2/ 4),  in.2 

Cb  =  drag  coefficient,  drag  force/(0.5 =  L"  Aref) 

A  eg  =  mass  offset  radial  distance  from  body  center,  in. 

Ci  =  rolling  moment  coefficient,  I/(  Arefd) 

C/j  =  induced  roll  moment  coefficient,  Equation  (1) 

C/p  =  roll  moment  damping  coefficient  derivative,  dCi/  d(pd/V),  per  radian 

C\S  =  fin  produced  roll  moment  coefficient  derivative,  dCil  d8 ,  per  radian 

C/v  =  normal  force  coefficient,  normal  force/(  Aref) 

Cjsra  =  normal  force  slope  coefficient,  dCN  /  da,  per  radian 

Cm  =  pitching  moment  coefficient,  pitching  moment/ ( Arefd) 

Cmo  =  asymmetry  pitching  moment  coefficient 

Cm  a  =  pitching  moment  slope  coefficient,  dCM  /  da ,  per  radian 

d  =  projectile  reference  diameter,  in. 

Ix  =  axial  (polar)  moment  of  inertia  about  the  body  spin  axis,  lb-in.2 

Iy  =  transverse  moment  of  inertia  about  an  axis  passing  through  the  body 

center  of  mass,  lb-in.2 
I  =  roll  moment,  lb-in. 

L  =  total  length  of  the  projectile,  in. 

M  =  Mach  number  of  the  projectile 

n  =  number  of  fins  in  a  fin  set 

rid  =  number  of  damaged  fins  in  a  fin  set 

p  =  spin  rate  of  projectile,  rad/s  (except  where  otherwise  noted  in  Hz) 

O' oo  =  dynamic  pressure,  (0.5/7  V2),  psi 

At  =  time  step  for  numerical  integration,  s 

V  =  projectile  velocity,  ft/s 


Greek  Symbols 


a  =  body  pitching  angle  of  attack,  deg 

=  d«  /  dt 

(  r~i — iA 

a  t  =  total  angle  of  attack  (yaw)  i  =  -y  a  +  ,  deg 


(5  -  body  side  slip  angle  of  attack,  deg 

SeCj  =  equivalent  fin  cant  angle  for  a  whole  fin  panel  for  the  partially  canted 

(chamfered)  fin,  deg 

8(i  -  design  value  of  the  fin  cant  angle,  deg 

S,iH  =  value  of  the  fin  cant  damage  angle,  deg 

S1,S2,S3  =  angles  of  the  damaged  area  of  the  fin  panel,  deg 
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Aj,A2,A3  =  fin  damage  parameters:  deg,  deg,  and  inch,  respectively 
y  -  fin  plane  roll  orientation  angle  with  respect  to  the  plane  of  angle  of 

attack,  deg 

r  =  initial  angle  of  mass  offset  location,  deg 

A  =  angle  of  the  Cmo  asymmetry  moment,  deg 

p  -  air  density,  slug/ft3 

0  =  roll  angle,  radian 

=  roll  (spin)  rate,  d^/dt,  rad/s  (except  otherwise  noted  in  Hz) 


a>n  =  first  natural  frequency  in  lateral  bending,  rad/s 

a>p  =  pitching  motion  frequency,  J(pCMf/q/\-efd  /  I y\ rad/s 
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